Dynamics and Interactions of Binaries and Neutron Stars in Globular Clusters

We model the dynamics of test binaries in isotropic, multi-mass models of galactic globular clusters. The evolution of binary orbits through the cluster potentials is modeled, including second order diffusion terms, and probabilities for close encounters with field stars are calculated. We carry out Monte Carlo simulations of the effects of the binary--single star encounters on the binary population and distribution in the cluster, and estimate the collision rate for different stellar populations in globular clusters with different structural parameters. Assuming a Salpeter IMF, for low concentration clusters the core encounter rate is dominated by turnoff mass main--sequence stars and medium mass white dwarfs. For high concentration, high density clusters the encounter probabilities are increasingly dominated by neutron stars and heavy white dwarfs. Hence we predict a smaller ratio of blue stragglers and cataclysmic variables to pulsars in high concentration clusters. The total number of millisecond pulsars, and the ratio of single to binary pulsars, is broadly consistent with the observed population, suggesting the binary--single star encounters contribute significantly to the pulsar formation rate in globular clusters, for the whole range of globular cluster types. The number of millisecond pulsars and the ratio of pulsars in different globular clusters is best explained by a total binary fraction comparable to that of the galaxy, and a modest number of primordial neutron stars in the globular clusters.Comment: 59 pages, uuencoded compressed postscript, including 18 figures. Astrophysical Journal Supplements, in pres

Close Approach during Hard Binary--Binary Scattering

We report on an extensive series of numerical experiments of binary--binary scattering, analysing the cross--section for close approach during interactions for a range of hard binary parameters of interest in globular cluster cores. We consider the implied rate for tidal interactions for different globular clusters and compare our results with previous, complementary estimates of stellar collision rates in globular clusters. We find that the collision rate for binary--binary encounters dominates in low density clusters if the binary fraction in the cluster is larger than $0.2$ for wide main--sequence binaries. In dense clusters binary--single interactions dominate the collision rate and the core binary fraction must be \ltorder 0.1 per decade in semi--major axis or too many collisions take place compared to observations. The rates are consistent if binaries with semi--major axes $\sim 100 AU$ are overabundant in low density clusters or if breakup and ejection substantially lowers the binary fraction in denser clusters. Given reasonable assumptions about fractions of binaries in the cores of low density clusters such as NGC~5053, we cannot account for all the observed blue stragglers by stellar collisions during binary encounters, suggesting a substantial fraction may be due to coalescence of tight primordial binaries.Comment: 13 pages including 13 ps figures. MNRAS in pres

The orbital statistics of stellar inspiral and relaxation near a massive black hole: characterizing gravitational wave sources

We study the orbital parameters distribution of stars that are scattered into nearly radial orbits and then spiral into a massive black hole (MBH) due to dissipation, in particular by emission of gravitational waves (GW). This is important for GW detection, e.g. by the Laser Interferometer Space Antenna (LISA). Signal identification requires knowledge of the waveforms, which depend on the orbital parameters. We use analytical and Monte Carlo methods to analyze the interplay between GW dissipation and scattering in the presence of a mass sink during the transition from the initial scattering-dominated phase to the final dissipation-dominated phase of the inspiral. Our main results are (1) Stars typically enter the GW-emitting phase with high eccentricities. (2) The GW event rate per galaxy is a few per Gyr for typical central stellar cusps, almost independently of the relaxation time or the MBH mass. (3) For intermediate mass black holes (IBHs) of ~a thousand solar masses such as may exist in dense stellar clusters, the orbits are very eccentric and the inspiral is rapid, so the sources are very short-lived.Comment: ApJ Accepte

Transition from adiabatic inspiral to plunge into a spinning black hole

A test particle of mass mu on a bound geodesic of a Kerr black hole of mass M >> mu will slowly inspiral as gravitational radiation extracts energy and angular momentum from its orbit. This inspiral can be considered adiabatic when the orbital period is much shorter than the timescale on which energy is radiated, and quasi-circular when the radial velocity is much less than the azimuthal velocity. Although the inspiral always remains adiabatic provided mu << M, the quasi-circular approximation breaks down as the particle approaches the innermost stable circular orbit (ISCO). In this paper, we relax the quasi-circular approximation and solve the radial equation of motion explicitly near the ISCO. We use the requirement that the test particle's 4-velocity remain properly normalized to calculate a new contribution to the difference between its energy and angular momentum. This difference determines how a black hole's spin changes following a test-particle merger, and can be extrapolated to help predict the mass and spin of the final black hole produced in finite-mass-ratio black-hole mergers. Our new contribution is particularly important for nearly maximally spinning black holes, as it can affect whether a merger produces a naked singularity.Comment: 9 pages, 6 figures, final version published in PRD with minor change

Resonant relaxation near a massive black hole: the stellar distribution and gravitational wave sources

Resonant relaxation (RR) of orbital angular momenta occurs near massive black holes (MBHs) where the stellar orbits are nearly Keplerian and so do not precess significantly. The resulting coherent torques efficiently change the magnitude of the angular momenta and rotate the orbital inclination in all directions. As a result, many of the tightly bound stars very near the MBH are rapidly destroyed by falling into the MBH on low-angular momentum orbits, while the orbits of the remaining stars are efficiently randomized. We solve numerically the Fokker-Planck equation in energy for the steady state distribution of a single mass population with a RR sink term. We find that the steady state current of stars, which sustains the accelerated drainage close to the MBH, can be up to ~10 times larger than that due to non-coherent 2-body relaxation alone. RR mostly affects tightly bound stars, and so it increases only moderately the total tidal disruption rate, which is dominated by stars originating from less bound orbits farther away. We show that the event rate of gravitational wave (GW) emission from inspiraling stars, originating much closer to the MBH, is dominated by RR dynamics. The GW event rate depends on the uncertain efficiency of RR. The efficiency indicated by the few available simulations implies rates ~10 times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting such events by future GW detectors, such as LISA. However, a higher, but still plausible RR efficiency can lead to the drainage of all tightly bound stars and strong suppression of GW events from inspiraling stars. We apply our results to the Galactic MBH, and show that the observed dynamical properties of stars there are consistent with RR.Comment: Accepted to ApJ; Minor revision

Gravitational radiation timescales for extreme mass ratio inspirals

The capture and inspiral of compact stellar masses into massive black holes is an important source of low-frequency gravitational waves (with frequencies of ~1-100mHz), such as those that might be detected by the planned Laser Interferometer Space Antenna (LISA). Simulations of stellar clusters designed to study this problem typically rely on simple treatments of the black hole encounter which neglect some important features of orbits around black holes, such as the minimum radii of stable, non-plunging orbits. Incorporating an accurate representation of the orbital dynamics near a black hole has been avoided due to the large computational overhead. This paper provides new, more accurate, expressions for the energy and angular momentum lost by a compact object during a parabolic encounter with a non-spinning black hole, and the subsequent inspiral lifetime. These results improve on the Keplerian expressions which are now commonly used and will allow efficient computational simulations to be performed that account for the relativistic nature of the spacetime around the central black hole in the system.Comment: 19 pages, 4 figures. Changed in response to referee's report. Accepted for publication in Astrophysical Journa

Gravitational waves from stars orbiting the Sagittarius A* black hole

One of the main astrophysical processes leading to strong emission of gravitational waves to be detected by the future space-borne interferometer LISA is the capture of a compact star by a black hole with a mass of a few million solar masses in the center of a galaxy. In previous studies, main sequence stars were thought not to contribute because they suffer from early tidal disruption. Here we show that, according to our simulations of the stellar dynamics of the Sgr A* cluster, there must be one to a few low-mass main sequence stars sufficiently bound to the central Galactic black hole to be conspicuous sources in LISA observations. The probability that a white dwarf may be detectable is lower than 0.5 and, in spite of mass segregation, detection of a captured neutron star or stellar black hole in the center of the Milky Way is highly unlikely.Comment: 5 pages, 3 figures, accepted for publication in ApJL, new version shortened to fit in 4 journal pages. Slightly longer version available at http://obswww.unige.ch/~freitag/papers/article_SgrA_long.ps.g